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ABSTRACT 

We present 0."1 resolution near-infrared integral field spectroscopy of Ha in a 
z=1.4781 star forming galaxy, Q2343-BM133. These observations were obtained 
with OSIRIS (OH Suppressing Infra-Red Imaging Spectrograph) using the W.M. 
Keck Observatory Laser Guide Star Adaptive Optics system. Ha emission is re- 
solved over a Of! 8 (6.8 kpc) x &!5 (4.3 kpc) region with a Of!l spatial resolution. We 
find a global flux of 4.2±0.6xl0~ 16 ergs s~ x cm -2 , and detect a spatially resolved 
velocity gradient of ~134 km s -1 across the galaxy and a global velocity disper- 
sion of 73±9 km s _1 . An upper limit of [Nll]/Ha< 0.12 is inferred, which implies 
that this galaxy is not dominated by an active galactic nucleus and has a metal- 
licity at or below 1/2 solar metallicity. We derive a star formation rate (SFR) 
of 47±6 M yr _1 , and a dereddened SFR of 66±9 M Q yr _1 . Two-dimensional 
kinematics for Q2343-BM133 fit well with an inclined-disk model, with which 
we estimate an enclosed mass of 4.3 xlO 9 M within 5.5 kpc. A possible merger 
scenario is also presented, and can not be fully ruled out. We derive a virial mass 
of 1.1 xlO 10 M for a disk geometry, using the observed velocity dispersion. We 
propose that Q2343-BM133 is currently at an early stage of disk formation at a 
look-back time of 9.3 Gyr. 

Subject headings: galaxies: evolution - galaxies: high-redshift - galaxies: kine- 
matics and dynamics - infrared: galaxies 
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Introduction 



Formation and evolution of galactic disks remains one of the most important unsolved 
aspects of extragalactic study. At redshifts up to z ~ 1, optical spectroscopy and Hub- 
wealth of data on disk galaxies (eg, 
However, at higher redshifts galaxy 



Voet et al. 


1996. 


Voet et al. 


1997. 


Lillv et al. 


1998) 



morphologies become difficult to study since their bolometric surface brightnesses fade as 
(1+z) 4 and rapidly shrink to angular scales under one arcsecond. Beyond a redshift of one, 
optical observations are only able to detect rest frame ultraviolet (UV) light and older stellar 
populations can be overlooked. Because of these observational limitations, spatial densities 
and properties of disk galaxies at higher redshift (z > 1) are still unknown, and we have 
little knowledge about the epoch of disk formation. 

With recent photometric and spectroscopic deep surveys our understanding of high 
redshift (1.4 < z < 3.5) star-forming galaxies has greatly increased, yielding properties of 
mass, stellar populations, gas and dust content, and star formation rates. Most successful 
among these, are color selection techniques including the search for galaxies with a strong 
Lyman Break between optical imaging filters. These Lyman Break Galaxies (LBGs) are 



(Steidel et al. 


1996; 


Madau et al. 


1996) 



early development of structures of star forming galaxies. To date, LBGs have primarily been 
studied in rest-frame UV light to investigate large-scale clustering, global star formation 
rates (SFR), and production of metals and outflows. HST observations of LBGs have shown 
structure down to 8 kpc in size, which is usually clumpy and irregular and is indicativ e of 



regions of active star formation (jConselice et aL||2003j ; lLotz et al.l 12006c lLaw et al.l 120071 ). 



There has been growing evidence which suggests that at least a s ubset of these LBG 
systems have disk-like morphology at rest-frame optical wavelengths (ILabbe et al.l 120031 ) 
and show large organized rota t ion which may indicate formation of an early galactic disk 
( Ivan Dokkum fc Stanfordll200ll . lErb et al. I l2003h . These studies have used NIR slit-spectroscopy 
of Ha emission under seeing-limited conditions. At these redshifts, the best seeing-limited 
observations probe roughly 5 kpc and yield only 2 or 3 resolution elements across the typi- 
cal galaxy. Many have argued that these kinematic measurements are insufficient to claim 
rotation since there are only three independent samples across the galaxies, and their curves 
could easily be duplicated with a variety of clumpy distributions. 

With the advent of integral field spectrographs (IFS), observations are now able to ex- 
tend this study to two-dimensional (2D) mapping of galaxy's kinematics and morphology. 
Recently, the IF S SINFONI (VLT) was used to conduct seeing-limited observations of 14 
galaxies at z~2 (IForster Schreiber et al.l 120061 ) . which showed four galaxies that exhibit ro- 
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tating disks with clumpy morphologies. However, the iForster Schreiber et al.l (120061 ) study 
also observed very high internal velocity dispersions (often > 100 km/sec) and concluded 
that their disks may be short lived and not precursors of modern disk galaxies. Again, with 
a seeing limited spatial resolution it is difficult to reliably interpret kinematics in objects 
that are roughly an arcsecond in diameter. A star forming galaxy has now also been imaged 
with SINFONI usin g a natural guide star adaptive optics (AO) system. This observation 
( jGenzel et al.l 120061 ) suggest that the z=2.38 galaxy has a massive rotating progenitor disk. 
However, their best fit model has residual velocities comparable to the original measure- 
ments, indicating that significant deviations from simple rotation are present. These IFS 
studies have targeted z > 2 galaxies, with a look-back time of more than 10 Gyr, which is 
greater than the expected formation age of the Milky Way disk. 

Using Keck Observatory's Laser Guide Star Adaptive Optics (LGS-AO) and the IFS 
of OSIRIS (OH Suppressing Infra-Red Imaging Spectrograph) we have obtained a spatially 
resolved (< O'.'l) spectrum of Ha emission from the z=1.4781 star forming galaxy Q2343- 
BM133. OSIRIS coupled with LGS-AO allows full two-dimensional mapping of Ha mor- 
phology and dynamics of the galaxy. We selected Q2343-BM133 as one of our targe t s from 
the rest-UV color selected catalog of Isteidel et~all J2003I . l2004h : kdelberger et al.l j2004h . 
Q2343-BM133 was chosen for its compactness, and because it is one of the bluest objects 
in the Steidel catalog with an R-K color of 2.1. In §2] we describe the observations and re- 
ductions. In $3] we present general Q2343-BM133 morphology and kinematics, and in $4] we 
discuss Ha kinematics, disk modeling, mass distribution, metallicity, and a possible merger 
scenario. In this paper we assume A-domin ated cosmology with Q M =0.27, Q v =0.73, and 
H =71 km sec -1 Mpc -1 (IBennett et al.ll2003l ). For our redshift of z=1.4781 this corresponds 
to ~ 8.53 kpc per arcsecond and a luminosity distance of 10.8 Gpc. 



Observations & Data Reduction 



Q2343-BM133 (hereafter: BM133) was the first successful detectio n of a high-redshif t 
star forming galaxy using the new integral field spectrograph OSIRIS (ILarkin et al.l 120061 ). 
We obtained data on BM133 with OSIRIS in conjunction with the Keck LGS-AO system 
( iWizinowich et al.ll2006l ; Ivan Dam et al.l 120061 ). as part of commissioning time on November 
23, 2005. The OSIRIS spectrograph utilizes a lenslet array, grating, and 2048x2048 Rockwell 
Hawaii II HgCdTe infrared detector to obtain simultaneous spectra over a rectangular field 
of view at a spectral resolution of R~3000. Reimaging optics in OSIRIS allow the user to 
change the spatial resolution on the lenslet array. We observed BM133 in the Hn3 (A cen = 
1.6348/im, AA = 0.0814 /mi) filter with a spatial sampling of O'.'l in both lenslet axes, which 
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provides a 4'.'8 x 6'.'4 field of view. 

The LGS-AO system produced an artifical guide star directly on the OSIRIS optical 
axis, with the tip-tilt sensor closed on an R = 16.2 mag star located 55'.'2 away from BM133. 
We chose an instrument position angle of 90° to maintain the science target and tip-tilt star 
within the AO and OSIRIS fields of view, respectively. To ensure acquisition and assess the 
quality of the night we observed the tip-tilt star with OSIRIS LGS-AO before moving to the 
galaxy, and measured the point spread function full-width at half-maximum (FWHM) of < 
O'.'l. This is an upper limit to the optical system performance since we have undersampled the 
PSF of the AO correction at this lenslet resolution. For the given magnitude and separation 
of the tip-tilt star, the Keck LGS-AO system predicts a Strehl ratio of 24%. For each 
dataset we took a 900 second frame on source with BM133 shifted 1" north of the field 
center, followed by a 5" eastern dither for a 900 second sky exposure, and then another on 
source 900 second exposure with BM133 shifted 1" south of the field center. We repeated 
this data set three times, yielding a total exposure of 90 minutes on BM133 with three 900 
second sky exposures to ensure good quality sky subtraction. 

Reductions were performed using the OSIRIS final data reduction pipeline (DRP) in 
combination with a custom IDL emission-line extraction routine. The DRP reduction steps 
used are similar to those performed for other NIR spectrograph instruments: correction of 
detector non-linearity, removal of bad pixels, removal of detector cross talk, sky subtraction, 
wavelength calibration, and mosaicing of images. Wavelength calibration is determined 
uniquely for each of the resulting 3000 spectra using arc line lamps. The critical and unique 
reduction steps for this data set are extraction of emission line spectra and assembly of the 
3D cube (x and y spatial locations of lenslets and A wavelength), in which the flux for each 
lenslet is the integral of the peak pixel values for each wavelength channel. 

3. Results 

The integrated Ha flux was collapsed from the final 90 minute mosaic frame using 
an average over a bandpass of AA = 0.014 /mi, which is presented in Figure [D The final 
mosaiced image has been smoothed spatially with a Gaussian kernel (FWHM 0'.'2) to achieve 
higher signal-to-noise for our kinematic analysis. The Ha flux emission of BM133 is resolved 
spatially and extends over a 0'.'8 (6.8 kpc) x 0."5 (4.3 kpc) region oriented from the northwest 
(NW) to the southeast (SE). We define the center of BM133 as the dynamical center as 
determined in §H and use this position for all figures. The morphology of BM133 is suggestive 
of a disk-like galaxy with Ha brightest in the NW region, and the Ha flux tailing off towards 
the SE region of the galaxy. This may indicate that the presence of dust extinction and/or 
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the distribution of Ha star-forming (SF) regions are not uniform throughout the galaxy, 
but may have a more clumpy distribution as seen in other high-z s tar forming galaxies 
dPapovich etaDbood : iForster Schreiber et al.lbood : benzel et al.ll2006h . 



Ha emission from BM133 has previously been d etected with NIRSPEC with a derived 
Ha flux of 3xl0 -16 ergs s _1 cm" " 2 ferb et al.ll2006al lbh. This was a seeing-limited slit spec- 
troscopy observation, which yielded a flux measurement with no 2D spatial information. 
Figure [2] is the spectrum of the summation of the entire spatial extent (1'.'3 x l'/3) of the 
unsmoothed final mosaic frame of BM133, where we find a total i ntegrated Ha flux of 



4.2±0.6xl0~ 16 ergs s" 1 cm" 2 . This is 40% more flux than found bv ferb et all (|2006al Jbh 



but we were not biased by slit orientation and are able to integrate the flux across the entire 
extent of the galaxy. Ha emission provides an instantaneous measurement of SF regardless 
of the global stellar population and history. We find an Ha luminosity of 5.9xl0 42 ergs s~ ] 
and u sing the Global Schmidt law (SFR (M yr" 1 ) 



Lf j r I 1 26xl0 41 ergs s" 



Kennicutt 



1998|) we find a global SFR of 47±6 Moyr" 1 . Erb et al J 12003 used optical colors and SED 
fitting to find an E(B-V)=0.115 for BM133. Applying this correction to the Ha luminosity 
results in a dereddened SFR dered = 66±9 M yr" 1 . We did not detect [Nil] (A rest = 0.6548 and 
0.6583 /iin) emission. However, using a square aperture of size 0'.'5 centered on the Ha flux 
location and the expected wavelengths of [Nil] emission, we were able to set an upper limit 
of [Nn]/Ha < 0.12. 

Fitting a Gaussian to the spatially integrated global spectrum of BM133 (Figure [2]) we 
find a systemic redshift of z= 1.4781. The Ha velocity width determined from the Gaussian 
width were corrected for the instrumental root-mean-square line width (43 km s _1 ) using 
atmospheric OH lines at nearby wavelengths, to yield a total velocity dispersi on (a v ) of 



73±9km s l . This is larger by 25 km s 1 compared to a v previously observed by lErb et al. 



f)2006al lbh with slit spectroscopy, and again our difference may be explained since OSIRIS 
samples the full Ha flux emission across the galaxy. This discrepancy might also be due to 
the difficulty to measure a v since NIRSPEC has a factor of 2 lower spectral resolution. Table 
1 summarizes these observational properties of BM133. 

To measure the kinematic structure of BM133 we fit a Gaussian for each lenslet with 
an integrated signal-to-noise (S/N) > 20, yielding a velocity offset relative to the systemic 
redshift, and a velocity dispersion per lenslet. To account for nearby atmospheric OH sky 
line residuals, we determined weights for our Gaussian fits from the standard deviation per 
wavelength channel on a sky region within the 3D cube. Figure [3] is a map of the velocity 
offsets at each lenslet, overlaid with a model spider diagram described in §4.11 We find a 
velocity range of ~134 km s -1 across the galaxy, with negative velocity offsets concentrated 
in the NW region and positive velocities within the SE region. Note that the level of shear 
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we see within BM133 would not have been detected in seeing-limited observations such 



as in the SINFONI dataset (IForster Schreiber et al.ll2006l ). The Ha velocity widths were 



also determined from the Gaussian fits at each lenslet location. The Ha widths were then 
corrected for the instrumental line widths. These final line widths include the effects of beam 
smearing from the wings of the PSF and from the range of velocities that are spanned by 
a lenslet in the smoothed data cube (approximately 0".2 lenslet -1 ). Therefore, we consider 
these widths as upper limits to the velocity dispersions present at each spatial location. The 
resulting map of intrinsic velocity dispersions (cr„) at each spatial location is presented in 
Figure HI and are further discussed in §4.11 



4. Analysis & Discussion 

4.1. Ha Kinematics 

The 2D velocity map for BM133 is indicative of a galaxy with a symmetrically rotating 
disk. This is seen in our velocity offsets with a velocity gradient across the major axis of 
BM133 and a steep, but continuous, velocity transition from the approaching and receding 
portions of the galaxy. We are able to map the line-of-sight radial velocities of BM133, 
but in order to know the true circular velocities of the galaxy we need to account for the 



orientation and projection effects. We have generated inclined-disk models (lBegemanlll987l ) 
and performed least squares fits to our observed velocity fields in order to find the appropriate 
disk model for BM133. The parameters in our model are rotation center of the galaxy (x ,y ), 
systemic velocity (V ), position angle of the major axis (</>), inclination angle (i), and circular 
velocity (V C (R)) as a function of radius from the center. We adopt a rotation curve profile 
that is commonly seen in local spiral galaxies, with a linearly rising circular velocity at lower 
radii and a plateau at larger radii. For V C (R) we use two parameters for the fit using a slope 
(m) and plateau radius (r p ) with the following equations, 



V c = mR, for R < i p 



V c = mr p , for R > r p (2) 

We modeled spider diagrams based on our seven parameters (x D , y Q , V D , i, 0, m, r p ) 
and convolved our models with the same Gaussian kernel (FWHM = Of! 2) as our observed 
data. We performed a grid search over all parameters to find the global minimum reduced 
chi-squared (x 2 ). This minimum occurs at x 2 = 0.72 with the parameters: z=54°, = 
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160°, T p = ft' 11 (0.9 kpc), and m = 65 km s _1 O'.'l" 1 . The corresponding model spider 
diagram is overlaid on our velocity field in Figure [3] with respect to the dynamical center of 
the galaxy. The center position (x ,y ), systemic velocity (V ), and position angle (0) were 
tightly constrained on the observed fields, with a significant degeneracy of the model found 
with the inferred circular velocity and inclination. This is expected since V C (R) and i are 
correlated by V c (R)sin(«), and t here is a known degeneracy between the two parameters that 
increases at larger inclinations (IBegeman! 119871 ) . A second local x 2 minimum occurs at i = 
86°. However, based on the ellipticity of the major and minor axes we rule out the second 
inclination (i = 86°) since the morphology is inconsistent with an edge-on galaxy. 

To compare with previous slit spectroscopy observations, in Figure [5] we have plotted 
the velocity profile for a projected slit across the major axis of the galaxy using our derived 
position angle of 0=160°. Weighted averages of the velocities were determined for each 
spatial position along the virtual slit using the flux S/N as weight for each velocity position. 
This illustrates the gradual velocity gradient across the SE portion of the galaxy, and the 
sharp transition and neg ative gradient acr os s the dynamical center to the NW region. As 



seen m previous studies ferb et al.l l2006al lbl: [ Forster Schreiber et al.ll2006l ) the Ha peak is 



not located at the dynamical center and can be explained by having the Ha emission arise 
in discrete SF regions spread through the disk, while the galactic center has a lower state of 
SF. 

Subtracting model velocities from our observed velocities results in residuals from 2 
- 35 km s" 1 , with a mean residual of 13 km s" 1 , which is far smaller than the observed 



residuals of ~170 km s _1 from the iGenzel et al.l (120061 ) rotating disk model for their z=2.38 
galaxy. Nearby spiral ga laxies show typical de viations of ~ 10 km s" 1 in their velocity profile 
caused by spiral ar ms (jSofue fc Rubinl 1200 ll ) and intrinsic dispersions within H// regions 
( jWeiner et al.ll2005l ). The AO PSF will further smooth the velocity profiles and increase 
the measured dispersions, since each spatial element of our data cube averages over the AO 
PSF per wavelength channel and the spectral features are broadened. However, this effect of 
beam smearing can not be fully determined without knowing the intrinsic distribution of Ha 
emission. We investigate beam smearing effects by simulating a galaxy's Ha flux distribution 
of both an exponential disk or two compact sources in an OSIRIS data cube, convolved with 
varying PSF performances of the AO system. For a Strehl ratio of 20% close to our projected 
AO performance, the typical beam smearing on velocity dispersions range from 18 to 33 km 
sec -1 for an exponential disk derived from our modeled inclined disk parameters. For two 
compact sources with differing intensities we find velocity dispersions from 5 to 55 km sec" 1 , 
where the peak dispersion occurs at the brightest knot. Both of these models, imply that a 
significant portion of the observed velocity dispersions are intrinsic to BM133, particularly 
within the NW region. 
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The mean internal velocity dispersion for a given spatial location in BM133 is 68 km 
s _1 with the highest a v near the flux center, ju st NW of the dynamica l cente r. Smaller a v 
values occur SW from the dynamical center. In iForster Schreiber et al.l (120061 ) much higher 
dispersions (often greater than 100 km s _1 ) are observed. Although their seeing limited 
observations likely suffer from significantly more beam smearing, this is unlikely sufficient 
to explain the differences between the two samples. They use the large ratio of velocity 
dispersion to circular velocity (<J v /v c ) to argue that the disks are temporary and will likely 
contract to form bulges as they dissipate random motion. Similarly, the single z=2.38 galaxy 
observed by lGenzel et al.l (120061 ) with AO also shows very high velocity dispersions wh i ch are 
too large to be explained through beam smearing. The best fit disk in iGenzel et al.l (120061 ) 
also leaves residuals of as much as 170 km s _1 , implying even stronger velocity components 
that are not from simple disk rotation. With our observed dispersions of 68 km sec -1 and 
mean residual velocities of 15 km s -1 of our simple disk model, we believe that BM 133 is a 
candidate for a true disk that may persist to the present day. In addition with the ability of 
LGS-AO we are able to resolve a compact knot region (~200 pc) in the NW of the galaxy, 
this may be indicative of high star formation region(s) either within a portion of the disk 
or the possiblity of an early-stage merger as we discuss in Section §4.41 We note that its 
redshift of 1.4781 places it at a look back time of 9.3 Gyr, comparable to the estimated age 
of the Milky Way Disk of 7.3±1.5 Gyr from studies of white dwarfs ([Hansen et al.l 12 002) 



and o f 7.5 — 10 Gyr from isochrone fitting of the thick-disk cluster NGC 6791 (jSandage et al. 
2003h . 



4.2. Mass Distribution 

Using our disk model and observed kinematics we are able to investigate the mass 
distribution of BM133. We calculate an enclosed dynamical mass (M enc i) of a highly flatten 
spheriod assuming Keplerian motion, 

2v 2 r 

M end = (3) 

For r, we choose the farthest distance from the dynamical center with a measured velocity 
offset at r = / .'64 (5.5 kpc). We use the circular velocity (v c = v p = 72 km s _1 ) determined 
from the inferred spider diagram, and the gravitional constant (G), to find a dynamical mass 
of M encl = 4.3 xlO 9 M . 

The enclosed mass is a lower limit on the total mass of BM133 since it underestimates the 
halo mass of the galaxy which includes stellar, gas, and dark matter mass components beyond 
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the radius (r) assumed. Following iForster Schreiber et al.l (120061 ). we are ab 



the halo mass (Mw» ) of BM133 using the spherical virialized collapse model ([Peebles! 11980 



e to estimate 



White fc FrenMll99lh . These models estimate the halo mass with the following equation, 



Mi 



halo 



V, 



2 V1.67 



H- 1 G- 1 (1 + 



-1.5 



(4) 



We find a total halo mass of M^aio = 3.4 xlO 10 M^for BM133. T his halo mass is roughly an 
order-of-magnitude lower than the ( Forster Schreiber et al. 2006h z > 2 samp le. In addition, 
we confirm our halo mass by deriving a similar estimate from lErb et al.l (120031 ) by recognizing 



that Qm/^b ~ 7, where VLb is the total baryonic matter from M ste u ar and M gas . Therefore, 
we place a limit on the total halo mass of BM133 to be 7 times larger than its stellar 
mass, which likely dominates within the Ha region. With our calculated enclosed mass we 
find Mhaio ^ 4.8 xlO 10 M , which is in reasonable agreement with our spherical virialized 
collapsed model estimate. The halo mass derived is a lower limit given that modern galaxies 
have much smaller ratios of baryonic to dark matter. 

The virial theorem provides an additional way of estimating the enclosed mass of a 
galaxy, assuming velocity dispersion is due to the collective motion of the gas and stars in a 
gravitionally bound system according to the following equation, 



G 



(5) 



The galaxy's mass distribution and velocity isotropy influences the constant factor (C), 
which usually ranges from 1 to 10. Most of these contributions are unknown, previous wor k 
adopted a uniform spherical dis tribution o f C=5 (jPettini et al.ll200ll ; lErb et al.ll2003l 120041 ) . 
For comparison, we follow the lErb et al.l (j2006al ) assumption of disk geometry and use C 
= 3.4. With an effective half light radius (r 1/2 ) for BM133's Ha flux of ~0"3 (2.56 kpc), 
and using the global yelocity di s persion we find a virial mass of M V i r = 1.1 xlO 10 M . This 
is comparable to the lErb et al.l (j2006al ) estimate of 1.8 xlO 10 M Q , assuming their observed 
a v and their assumed galaxy radius (~8 kpc). Our virial mass estimate is a factor of 2.6 
greater than our estimated enclosed mass. Several factors likely explain this relatively large 
disagreement. In particular, by using (7=3.4 we are assuming an inclination of 45° instead of 
our model's value of i=54°. More importantly, the dispersions are likely significantly higher 
than the model assumes due to the irregular distribution of Ha flux being a biased tracer of 
the random motion. In addition, we can not rule out a merger scenario where two knot-like 
structures are infalling, and therefore induce high SF in each merger component. A possible 
merger scenario is further discussed in §4.41 
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Of all the mass estimates we favor the enclosed mass from the inclined-disk model since 
it includes inclination effects and uses intrinsic circular velocities. However, we recognize that 
this enclosed mass is substantially lower than other mass estimates of star forming galaxies at 
this redshift range. Our mass estimates are probing the inner kiloparsecs of the galaxy, and 
BM133 is likely to be imbedded i n a mu ch larger, extended structure beyond our observed 
radius. We note that lErb et al.l (j2006al ) find a gas fraction mass of ~85% giving BM133 
one of the highest star formation rates per unit stellar mass in the NIRSPEC sample. This 
large gas fraction will likely lead to further contraction which will steepen its velocity profile 
and produce larger velocities overall. Therefore the slope of the observed velocity profile is a 
lower limit to the eventual modern value. The kinematics observed still fit within the range 
of rotation curves found in nearby spiral galaxies with comparable enclosed massses out to 



these radii (ie, NGC 598 : ISofue et al.lll999l : ISofue fc Rubinll200l[ ). However, we stipulate 



that BM133 may continue to grow into a more massive disk system given its observed youth 
and gas richness, where the majority of the growth will occur at larger radii. 

The enclosed, virial, and halo masses inferred here are model dependent and are ex- 
trapolations from the observed Ha distribution and kinematics which are concentrated near 
the center of the galaxy. We investigate th ese quantities in part to compare with prior 
mass estimates of BM133 (|Erb et al.ll2006al ) and mass estimates of similar objects in the 
Forster Schreiber et al.l (120061 ) survey. A significant difficulty in such a comparison is that 
the IFS LGS-AO observations focus on the inner several kpc, while the seeing-limited cam- 
paigns extend to larger radii, albeit with poorer resolution. Therefore the mass estimates 
will likely represent a lower limit since we do not probe the entire distribution of even the 
baryon matter of the galaxy. 



4.3. Metallicity 



We did not detect optical rest-frame [Nil] emission, however we were able to place an up- 
per limit to [Nil] /Ha < 0.12. With suc h a low [Nil] /Ha rat io we can rule out that BM133 is 
dominated by an active galacic nucleus ([Baldwin et al.lll98ll ) . Using the r elationship between 



oxygen abundance ( O/H) and [Nil] /Ha p reviously used in LBG studies (IPettini et al.ll2001 



Shapley et al.ll200ll ; IPettini fc Pagelll2004j ) we can estimate the H// region abundances. With 
this relationship we find an oxygen abundance upper limit of 12 + logfO/ H) < 8.38. The 
assum ed oxygen abundance for solar metallicities is 12 + log(0/H) = 8.65 (IPettini fc Page! 
2004 ). therefore we infer a metallicity less than 1/2 solar for BM133. Our [Nil] /Ha limit is 



consi s tent w i th other abundance studies for z>2 LBG samples (jShapley et al.ll200ll ; lErb et al. 
20031 . l2006d ; iForster Schreiber et al.l 120061 ) . Our inferred (O/H) abundance limit is in rea- 



sonable agreement with lErb et al.l (j2006d j relantionship between (O/H) abundance versus 
stellar mass, assuming their stellar mass of 5xl0 9 M from SED fitting. 



4.4. Possible Merger Scenario 

Although a smooth disk-like velocity structure can explain the observed kinematics, we 
note that with the LGS-AO resolution, much of the Ha flux is found in a single compact 
region (~200 pc) to the NW of the kinematic center. It is difficult to determine if this is a 
region of extremely high star formation within a disk, or a separate merging galaxy. In either 
case, the intensity of star formation in such a small region is extreme with an observed SFR 
> 16 M yr _1 , and the region shows the highest internal dispersion within all of BM133. If 
this is a star forming region within an otherwise stable disk, then it represents an extreme 
mode of star formation. But in this gas rich system it might be possible to temporarily 
support such run-away star formation. It seems reasonable that such a region would have 
significantly higher internal turbulence compared to modern star forming regions helping to 
explain the high intrinsic dispersion. 

In order to investigate this scenario we measure velocity dispersions (a v ) for two 0'.'3 x 
0'.'3 regions NW and SE of the galaxy. Assuming that each of these components reflect the 
internal dispersion of gas and stars and have a half-light radius of 11/2 = O'.'l (0.9 kpc), we 
find a virial mass for each component of M. V i r ,NW — 6.5xl0 9 M Q and M vil . t sE — 2.3xl0 9 M . 
If both of these components are infalling towards each other, we can estimate the expected 
observed velocity (v^sm i) by considering Keplerian motion with the given separation (3.81 
kpc) and virial masses. We find that for such a merger system the observed velocity would 
have an upper limit of 80 km s -1 . This velocity is quite similar to our observed plateau 
velocity of v p = 72 km s _1 , therefore we can not rule out the possibility of a merger scenario. 
However, if this was a merger it would also be fortuitous to have tidal stripping between 
the two components that would produce such a smooth velocity gradient with a dispersion 
highest near one of the flux centers as we observe in BM133. So, although we cannot fully 
rule out an early-stage merger system, we believe that it is a less probable scenario given 
our measured two-dimensional velocity offsets and dispersions, and our well-fit inclined-disk 
model. 
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5. Summary 

We have presented first extragalactic results from the new integral field spectrograph 
OSIRIS with the Keck LGS AO system by measuring Ha emission from a z=1.4781 star 
forming galaxy. We spatially resolve Ho; over a 6.8 x 4.3 kpc region and find a global flux 
of 4.2±0.6xl0~ 16 ergs s" 1 . We find a high dereddened SFR of 66±9 M yr _1 . Based on 
[Nil] /Ha ratio limit we infer a 1/2 solar or below metallicity, which is consistent with other 
high-z LBG samples. We fit an inclined-disk model to the 2D kinematics of BM133, and 
find low residuals with a disk having a position angle (0) of 160°, inclination angle (i) of 54°, 
and plateau velocity of 72 km s -1 achieved within a radius of ~ O'.'ll (0.9 kpc). We estimate 
an enclosed mass of 4.3xl0 9 M and virial mass of 1.1 xlO 10 M , which reconfirms results 
by lErb et al. (|2006al Jbh. We cannot conclusively eliminate a merger scenario, but we favor a 



true disk in an early stage of star formation. 

An IFS combined with LGS-AO offers a powerful new tool for studying morphology and 
kinematics of z > 1 galaxies. The capability to achieve diffraction-limited NIR observations, 
comparable to optical HST observations, offers a superb and necessary capability for studying 
kinematics and evolution of high-z galaxies. Although results presented in this paper are for 
only one target, the analysis demonstrated here is being applied to a larger on-going study 
of star forming galaxies at these redshifts. The capability of near-IR IFS instruments such 
as OSIRIS and SINFONI to provide 2D mapping of optical rest-frame spectral features of 
high-z galaxies has already demonstrated significant advantages for investigating kinematics 
and mass distribution of these objects. BM133 has offered a stepping stone for future high- 
z extragalactic kinematic studies, and offers one of the best candidates for a precursor of 
present-day spiral galaxies. 
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Table 1. Q2343-BM133 Observational Properties 



RA 


DEC z Ha 




Exposure 


Date 


[Nil] /Ha 


Av 




(J2000) 


(J2000) 


(mags) 


(seconds) 


Observed (ergs s - 1 cm -2 ) 




(km s" 1 ) 


(km s _1 ) 


23:46:16.18 


12:48:09.31 1.4781 


20.5 


6x900 


November 2005 4.2±0.6xl0~ 16 


< 0.12 


134 


73±9 



t Observed bv lErb et all J2006al fbh 



Table 2. Q2343-BM133 Derived Properties 



SFR 

(Meyr- 1 ) 


SFR^ered 

(M a yr- 1 ) 


ri/2 
(kpc) 


(km s" 1 ) 


(kpc) 


Position Angle 

w 


Inclination 

(0 


M enc ; 
(M Q ) 


(M s ) 


M haio 
(Mq ) 


47±6 


66±9 


2.6 


72 


0.9 


160 


54 


4.3xl0 9 


l.lxlO 10 


3.4xl0 10 



tv p is used for v c in the enclosed and halo mass estimates 
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Fig. 1.— A Gaussian smoothed (FWHM=0'/2) mosaiced image of Q2343-BM133 with a total 
exposure of 90 minutes collapsed around Ha (AA = 1.4 nm) with a spatial size of 2 / .'0x2'. / 
and position angle of 0°. 
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Fig. 2. — Total Ha emission spectrum of Q2343-BM133, collapsed over a spatial size of 
l''3xl'/3 encompassing the entire extent of the galaxy. A Gaussian is fitted over the Ha emis- 
sion to yield our global velocity dispersion (cr v ) with intrumental correction performed. 
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ARA (") 

Fig. 3. — Radial Velocity: Two dimensional Ha kinematics of Q2343-BM133 showing the 
spatial distribution of velocity (km s" 1 ) relative to the measured systemic velocity. This was 
obtained from Gaussian fits to each spatial element (lenslet) after smoothing spatially with 
a 2D Gaussian of FWHM=2 (or &!2). The velocities are color-coded matching the labeled 
color-bar. Overlaid is the spider diagram modeled in £JH with each contour representing 10 
km s _1 . 
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Fig. 4. — Velocity Dispersion: Two dimensional Ha kinematics of Q2343-BM133 showing 
spatial distribution of velocity dispersion(km s _1 ) after removal of instrumental widths. 
This was obtained from Gaussian fits to the data cubes after smoothing spatially with a 
two-dimensional Gaussian of FWHM=2 (or 0'.'2). Velocities are color-coded according to 
the labeled color bar. Overlaid is the spider diagram modeled in £JH using the radial velocity 
offsets (see Figure [3]). Note we see higher velocity dispersions in the NW region of the galaxy, 
nearest the dynamical center, and lower velocity dispersions nearest the plateau velocities, 
indicated by the overlaid spider diagram with each contour representing 10 km s _1 . 
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Fig. 5. — Velocity (km s _1 ) versus the projected axis (") for the 160° position angle de- 
termined from the spider diagram modeling in §4.11 This plot represents a slit along the 
major axis of the galaxy. Projected velocty at each distance was determined from a weighted 
average of the velocities, using the flux signal-to-noise as weight. The center position (0.0) 
is defined from the dynamical center of the inclined disk model in §H Errors get larger 
with greater distance from the center since there were less spatial elements to weight in the 
average. 



